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Abstract 

This  study  presents  a  computational  fluid  dynamics  (CFD)  model  for  modelling  gas  evolution  and  current  distribution  in  a  direct  methanol  fuel 
cell  (DMFC).  The  improved  two-phase  model  includes  a  new  sub-model  for  estimating  the  interface  mass  transfer  without  empirical  correlations. 
Simulation  results  in  a  horizontal  channel  of  the  DMFC  agree  with  typical  trends  reported  in  the  literature  for  bubbly  flows.  The  increase  in  inlet 
flow  rate  is  found  to  lead  to  a  decrease  in  the  gas  content  in  the  outlet  of  the  anode  channels.  A  case  study  illustrates  applications  of  the  CFD  model 
for  modelling  gas  evolution  and  current  distribution  in  a  DMFC  with  a  parallel  flow-field  design.  Simulation  results  with  a  improved  two-phase 
model  provide  an  explanation  of  experimental  observations  of  a  transparent  DMFC  with  parallel  channels.  An  improved  three-dimensional  CFD 
model  includes  all  relevant  phenomena  and  is  valuable  for  gas  management  in  a  DMFC  design. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Understanding  a  two-phase  phenomenon  is  important 
for  designing  a  high-performance  direct  methanol  fuel  cell 
(DMFC).  Critical  issues  for  improving  the  DMFC  performance 
are  methanol  crossover,  gas  management  on  the  anode  side, 
and  water  management  on  the  cathode  side.  Gas  manage¬ 
ment  is  especially  important  in  a  DMFC  design  since  the 
methanol  electrochemical  oxidation  produces  carbon  dioxide 
flux  on  the  anode  side.  Removing  CO2  bubbles  is  required  to 
avoid  blocking  anode  channels  that  may  lead  to  limited  mass 
transport. 

A  number  of  physicochemical  phenomena  take  place  in  a 
DMFC,  including  momentum  and  mass  transfer,  electrochem¬ 
ical  reactions,  and  gas-liquid  flow  in  the  anode  and  cathode 
channels.  All  these  processes  are  coupled  and  result  in  a  need  for 
optimum  cell  design  and  optimum  operating  conditions.  Thus, 
good  understanding  of  these  complex,  interacting  phenomena  is 
essential  in  fuel  cell  design. 
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Researches  devoted  to  exploring  mass  transfer  and  electro¬ 
chemical  reactions  in  DMFCs  are  numerous  [1-11].  Scott  et 
al.  [1]  examined  the  feasibility  of  using  stainless-steel  mesh 
materials  as  flow  beds.  They  reported  electrochemical  perfor¬ 
mance  and  gas  management  characteristics  with  flow-beds  from 
flow  visualization  studies  on  the  anode  side.  Geiger  et  al.  [2] 
used  neutron  radiography  to  investigate  of  gas-evolution  pat¬ 
terns  in  anode  flow-fields.  It  was  found  that  gas  accumulates 
to  a  large  extent  at  the  inner  section  of  spiral  channels  and 
thereby  blocks  a  considerable  part  of  the  active  area.  They  noted 
that  a  spiral  type  of  flow-field  is  not  appropriate  for  the  anode. 
Argyropoulos  et  al.  [3]  used  acrylic  cells  to  investigate  visually 
gas  evolution  in  an  operating  DMFC.  They  studied  the  effect 
of  operating  parameters  and  flow-bed  design  on  gas  manage¬ 
ment.  It  was  concluded  that  an  increase  in  inlet  flow  rate  is 
beneficial  to  gas  removal.  Using  a  5  cm2  transparent  cell,  Lu 
and  Wang  [4]  investigated  the  effects  of  backing  pore  structure 
and  wettability  on  cell  polarization  characteristics  and  two- 
phase  flow  dynamics.  They  found  that  an  anode  backing  layer 
of  uniform  pore  size  and  high  hydrophilicity  is  preferred  for 
gas  management  in  the  anode.  Tuber  et  al.  [5]  compared  the 
performance  of  PEMFCs  and  DMFCs  with  serpentine,  parallel 
and  new  fractal  flow-fields.  The  results  indicated  that  serpen- 
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Nomenclature 

C  mass  fraction  (kg  kg- 1 ) 

D  diffusion  (m2  s- 1 ) 

F  Faraday  constant  (C  mol- 1 ) 

g  acceleration  (ms-2) 

H  membrane  thickness  (m) 

h  channel  height  (m) 

I  current  density  (A  m-2) 

Iq  exchange  current  density  (A  m-2) 

Ie  ionic  current  density  (A  m-2) 

7l  capillary-diffusional  flux  of  the  liquid  phase 
(kgm-2  s-1) 

k  permeability  of  porous  material  (m2) 

K  distribution  of  the  components 

L  molar  flow  rate  (mol  s - 1 ) 

M  molecular  weight  (kg  mol- 1 ) 

N  mass  flux  (kgm-2  s-1) 

n  number  of  electrons 

n&  electro-osmotic  drag  coefficient 

p  pressure  (Pa) 

s  stoichiometric  coefficient 

S  area  (m2) 

T  temperature  (K) 

u  velocity  (ms-1) 

Uq2  thermodynamic  equilibrium  potentials  of  oxygen 
reduction  (V) 

jyMeOH  thermodynamic  equilibrium  potentials  of 
methanol  oxidation  (V) 
v  water  velocity  (ms-1) 

Vceii  cell  voltage  (V) 

Vanode  volume  of  anode  channels  (m3) 
v  molar  fraction  in  liquid  phase  (mol  mol- 1 ) ;  coor¬ 

dinate,  (m) 

y  molar  fraction  in  gas  phase  (mol  mol-1);  coordi¬ 

nate,  (m) 

z  coordinate  (m) 

Greek  symbols 

rQ  source  of  mass  in  gas  phase  (kg  m-3  s-1) 
a  a  charge-transfer  coefficient  of  the  anode 

ac  charge-transfer  coefficient  of  the  cathode 

£  porosity  (m3  m-3) 

sq  gas  content  (m3  m-3) 

y  local  fractional  vaporization;  kinetic  factor 

rj  overpotential  (V) 

ip  potential  (V) 

p  viscosity  (Pa  s) 

yc  advection  correction  factor 

p  density  (kgm-3) 

o  conductivity  (m) 

\[f  coefficient 

Subscripts 
i  component 

in  inlet 


out 

outlet 

L 

liquid 

G 

gas 

A 

anode 

C 

cathode 

eff 

effective 

mix 

mixture 

DL 

diffusion  layer 

ref 

reference  value 

t 

total 

s 

interface;  solid 

m 

membrane 

Superscript 

k  component  (MeOH,  CO2,  H2O,  O2) 


tine  flow-fields  give  both  the  highest  and  the  most  stable  perfor¬ 
mance. 

Since  it  is  very  difficult  to  measure  directly  concentration 
and  gas  content  profiles  in  anode  channels,  modelling  is  used 
to  study  a  DMFC.  Argyropoulos  et  al.  [6]  developed  a  model 
to  predict  the  local  pressure  and  chemical  composition  in  the 
anode  and  cathode  sides  of  a  liquid-feed  DMFC.  Birgersson  et 
al.  [7]  presented  an  isothermal,  two-dimensional,  liquid  phase 
model  for  the  conservation  of  mass,  momentum  and  species 
in  the  anode  channel  and  porous  media  of  a  DMFC.  The  data 
demonstrated  the  relative  importance  of  mass  transfer  resistance 
in  both  the  flow  channel  and  the  adjacent  porous  backing.  Mugia 
et  al.  [8]  derived  a  multi-component,  steady-state,  model  based 
on  phenomenological  transport  equations  for  the  catalyst  layer, 
diffusion  layer  and  membrane.  To  understand  the  role  of  model 
parameters,  they  performed  a  parametric  study  of  the  model 
together  with  experimental  validation.  A  comprehensive,  two- 
dimensional  model  of  two-phase  flow  with  multi-component 
transport  and  electrochemical  reactions  was  reported  by  Wang 
and  Wang  [9]  for  a  liquid-feed  DMFC,  including  electrodes, 
channels  and  PEM  separator.  Kulikovsky  et  al.  [10]  constructed 
a  two-dimensional  model  for  a  gas-fed  DMFC  with  a  new 
type  of  current-collector.  Schultz  and  Sundmacher  [11]  devel¬ 
oped  a  one-dimensional,  dynamic  model  of  a  DMFC  based  on 
Maxwell-Stefan  mass  transport  equations  and  a  Flory-Huggins 
activity  model. 

Optimum  flow  field  design  is  important  for  improving  flow 
patterns  and  gas  evolution  in  anode  channels.  CFD  simula¬ 
tion  is  widely  used  for  PEMFC  flow-field  design  [12-15]. 
Numerical  modelling  provides  a  better  understanding  of  the 
main  phenomena  that  govern  fuel  cell  performance.  Three- 
dimensional  modelling  is  important  to  capture  performance- 
limiting  effects  such  as  mass  transfer  and  gas  evolution.  The 
conventional  CFD-based  model  of  DMFC  requires  experimen¬ 
tal  correlations  for  closure  of  multiphase  model  equations 
prior  to  numerical  solution.  Empirical  correlations  limit  appli¬ 
cation  of  conventional  sub-models  for  gas-liquid  flow  in  a 
DMFC. 
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Our  previous  study  [16]  focused  on  a  model  development 
for  gas  management  in  the  anode  flow-field  of  a  DMFC  under 
the  assumption  of  uniform  current  density.  The  objective  of  this 
research  is  to  modify  the  three-dimensional  DMFC  model  to  pre¬ 
dict  the  cell  performance,  concentration  and  potential  profiles 
taking  into  account  two-phase  phenomena  in  anode  channels. 
The  improved  DMFC  model  includes  a  new  sub-model  for 
interface  mass  transfer  developed  without  using  empirical  cor¬ 
relations. 

2.  Model  formulation 

The  model  subdivides  the  fuel  cell  into  seven  regions,  namely: 
anode  flow  channel,  anode  diffusion  layer,  anode  electrocatalyst, 
membrane,  cathode  electrocatalyst,  cathode  diffusion  layer  and 
cathode  flow  channel.  Each  sub-model  is  described  below. 

2.7.  Two-phase  model  for  channels 

The  two-phase  model  is  widely  used  for  modelling  momen¬ 
tum  and  mass  transfer  in  gas-liquid  flows  [9,17-19].  Using  a 
flow  visualization  technique,  Yang  et  al.  [20]  developed  the  flow 
regime  maps  for  a  channel  with  a  gas  permeable  sidewall.  Bewer 
et  al.  [21]  reported  that  gas-liquid  flow  is  homogeneous  during 
a  homogeneous  bubble  discharge.  Triplett  et  al.  [19]  compared 
one-dimensional  simulation  results  with  experimental  data  and 
concluded  that  the  homogeneous  model  held  true  for  bubbly  flow 
in  channels.  As  mentioned  above,  Wang  and  Wang  [9]  applied 
a  homogeneous  two-phase  model  for  the  anode  channel  of  a 
DMFC.  Based  on  flow  visualization  data  and  real  application  of 
one-dimensional  and  two-dimensional  models,  this  study  con¬ 
siders  three-dimensional  gas-liquid  flow  in  anode  and  cathode 
channels  with  the  following  assumptions: 

•  the  fuel  cell  is  operating  at  a  steady-state, 

•  the  two-phase  flow  is  isothermal,  incompressible  and  homo¬ 
geneous, 

•  single-phase  flow  in  the  cathode  channel  is  isothermal  without 
evaporation  and  condensation. 

The  model  equations  in  Table  1  describe  the  distribution  of 
velocity,  the  gas  volume  fraction  and  the  mixture  concentrations 
in  the  anode  and  cathode  channels. 


Table  1 

Model  equations  for  anode  and  cathode  channels 


Governing  equations 

Mathematical  expression 

Continuity  equationa 

V • (pu)  =  0 

(1) 

Momentum  equation3 

V  •  (puu)  —  -  V p  +  V  •  T  +  pg 

(2) 

Stress  tensor3 

( dui  duj  2o  dun  \ 

(3) 

Continuity  equation 

v  •  ( SQpQUG )  =  Tq 

(4) 

for  gas  phase3 

Species  conservation13 

V  ■  0 ouCk )  =  V  •  ((1  -  sa)pLDkL  effVCkL 

+ 

SGPG£)G,effVCG) 

(5) 

3  Sokolichin  and  Eigenberger  [18]. 
b  Wang  and  Wang  [9]. 


Table  2 


Model  equations  for  diffusion  layer 


Governing  equations 

Mathematical  expression 

Continuity  equation3 

V • (pu)  =  0 

(10) 

Momentum  conservation3 

m  =  — (yp  +  pkg) 

(11) 

Ml 

Species  conservation 

V  •  (y,puCk)  =  V  •  (PL<effVC!  + 

pG^,effVCo)  -  V  '  «C£  -  CgWl) 

(12) 

Ohm’s  lawb 

VKeffV^s)  —  ^ 

(13) 

a  Wang  and  Wang  [9]. 
b  Kulikovsky  et  al.  [10]. 


For  single-phase  gas  flow  in  cathode  channels,  the  gas  vol¬ 
ume  fraction  is  sq  =  1  •  The  homogeneous  model  for  channels 
assumes  that  the  phases  move  with  the  same  velocity.  The  dif¬ 
ficulty  in  modelling  includes  multi-component  mass  transfer 
taking  place  across  the  interface.  The  mixture  variables  and  prop¬ 


erties  are  given  by: 

density  (p)  =  sqPg  +  (1  -  £g)Pl  (6) 

concentration  (pC)  =  CqsgPG  +  Cl(1  -  £g)Pl  (7) 

velocity  (pu)  =  uqEqPg  +  ml(1  -  £g)Pl  (8) 

viscosity  (/zeff)  =  /xG£g  +  j^L(l  “  SG )  (9) 

2.2.  Two-phase  model  for  diffusion  layer 


According  to  Wang  and  Wang  [9],  a  two-phase  model  for  the 
porous  diffusion  layer  is  given  by  the  model  equations  listed  in 
Table  2. 

Details  of  model  variables  and  mixture  properties  are 
described  by  Wang  and  Wang  [9]. 

2.3.  Electrokinetics 

Fluid  flow  and  mass  transfer  processes  occur  in  a  DMFC 
in  conjunction  with  electrochemical  reactions.  These  processes 
have  a  significant  impact  on  gas  management  and  fuel  cell  per¬ 
formance.  The  current  distribution  model  assumes  that: 

•  the  catalyst  layer  is  sufficiently  thin  to  treat  it  as  an  interface, 

•  isotropic  porous  media  exists  in  the  diffusion  layer  and  mem¬ 
brane, 

•  negligible  contact  resistance, 

•  fully  hydrated  membrane, 

•  Tafel  kinetics  for  irreversible  reactions. 

The  following  electrochemical  reactions  take  place  at  the 
anode  and  cathode  catalyst  layers,  respectively. 

CH3OH  +  H20  -*  C02  +  6H+  +  6e~  (14) 

02  +  4H+  +  4e“  —*■  2H20  (15) 

The  current  distribution  model  equations  account  for  elec¬ 
trochemical  reactions,  ohmic  losses  and  mass  diffusion  in  the 
diffusion  layer  and  catalyst  layer  regions  [10]  (Tables  3  and  4). 
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Table  3 

Model  equations  for  catalyst  layer 

Variable  Diffusion  layer/membrane  interface 


from  interface  mass  transfer  of  carbon  dioxide  in  the  anode  chan¬ 
nels.  Wang  and  Wang  [9]  used  a  mass  transfer  equation  with  an 
empirical  coefficient  for  estimating  the  source  of  the  mass 


Solid  potential 

Membrane 

potential 

Speciesa 

The  anode 
reaction  rate 

The  cathode 
reaction  rate 


OVeff ' 


'Vr^eff ' 


lz 

Q(Ps 

dz 


1 

-I 


(16) 

(17) 

(18) 

(19) 

(20) 


a  /eff  =  I  +  /p. 

Model  Eqs.  (16)-(22)  with  auxiliary  equations  and  bound¬ 
ary  conditions  are  taken  from  the  original  model  developed  by 
Kulikovsky  et  al.  [10]  and  Wang  and  Wang  [9].  Physicochemi¬ 
cal  properties  and  molecular  transfer  coefficients  are  calculated 
from  equations  cited  by  Wang  and  Wang  [9]. 

2.4.  Boundary  conditions 


„  A^,gQj,s  +  pg^gACg 

Eg  =  - 7 -  (24) 

h 

where  Nt&  is  the  total  mass  flux  transferred  from  liquid  to  gas 
phase;  /3q  the  mass  transfer  coefficient  in  the  gas  phase;  ACg 
the  driving  force  of  mass  transfer;  Cg,s  is  the  concentration  at 
the  interface.  Estimation  of  the  interface  mass  flux  with  the  non¬ 
equilibrium  model  requires  additional  sub-models  or  empirical 
correlations  for  mass  transfer  coefficients  in  the  gas  and  liquid 
phases.  The  mass  transfer  coefficient,  in  turn,  is  a  complicated 
function  of  hydrodynamics. 

As  shown  by  Argyropoulos  et  al.  [6]  and  Sundmacher  and 
Scott  [22,23],  the  equilibrium  flash  equation  is  suitable  for  com¬ 
puting  the  gas  content  in  the  anode  compartment  treated  as  a 
continuous  stirred  reactor.  To  define  the  interface  flux  from  liq¬ 
uid  to  gas  phase,  it  is  necessary  to  consider  the  equilibrium 
condition  in  the  multi-component  gas-liquid  flow  in  the  anode 
channels.  The  multi-component  mixture  includes  carbon  diox¬ 
ide,  methanol  and  water.  For  control  volume  8  V,  we  define 
fractional  vaporization  as  follows: 


Set  conditions  are  required  at  all  boundaries  of  the  computa¬ 
tional  domains.  At  the  inlet  of  both  the  anode  and  the  cathode 
flow  channels,  the  boundary  values  are  prescribed  from  the  sto¬ 
ichiometric  flow  rate  and  mass  fractions.  On  all  walls,  a  no-slip 
boundary  condition  is  applied  for  the  momentum  equations  and  a 
no-flux  boundary  condition  for  the  components.  Boundary  con¬ 
ditions  at  internal  interfaces  for  the  current  model  are  described 
in  detail  by  Mugia  et  al.  [8]  and  Kulikovsky  et  al.  [10]. 


where  8G  is  the  molar  flow  rate  transferred  from  liquid  to  gas 
phase;  L  the  molar  flow  rate  of  liquid  in  a  channel  with  volume 
8  V;  8 S  is  the  area  and  8 S  =  8 V/h. 

Using  the  fractional  vaporization  Eq.  (25),  we  suggest  the 
following  equation  for  estimating  the  source  of  mass 


2.5.  Gas  content  in  anode  channels 

The  basic  two-phase  model  is  given  by  Eqs.  (l)-(5).  The 
above  model  equations  are  coupled  closely,  so  the  whole  set 
must  be  solved  simultaneously  and  iteratively.  The  presence  of 
bubbles  in  the  gas-liquid  flow  is  reflected  by  the  gas  content  sq 
and  the  source  term  rG  that  accounts  for  interface  mass  transfer. 

In  accordance  with  the  two-phase  model,  the  gas  content  is 
found  from  the  continuity  equation  (4).  The  source  of  mass  in 
Eq.  (4)  by  definition  is: 

8G 

rG  =  MG  —  (23) 

Carbon  dioxide  is  produced  in  the  anode  channels  by  the 
liquid-phase  electrochemical  reaction  (14).  Gas  evolution  results 

Table  4 

Model  equations  for  membrane 

Governing  equations  Mathematical  expression 

Ohm’s  lawb  V(crmjeff,  V<pm)  =  0  (21) 

„  _»  k  n&M  7e 

Momentum  conservation*1  u  = - (V/?  +  pkg )  H -  (22) 

Ml  P  F 

a  Wang  and  Wang  [9]. 
b  Kulikovsky  et  al.  [10]. 


ra  =  Mar  (t  +  (26) 

where  Mq  is  the  molecular  weight  of  the  gas  phase;  \j/  the  coef¬ 
ficient,  \jr  =  L/ 8  V.  Local  fractional  vaporization  y  is  found  from 
solving  the  equilibrium  flash  equation: 

y^(X;  —  l)C'(Mmix/M,) 

E  (Kj  —  l)y  +  1  } 

where  K[  is  the  distribution  of  each  component  between  the 
vapour  and  liquid  phases,  Kt  =  ydxi;  y  is  the  local  fractional 
vaporization.  The  distribution  of  the  mixture  concentration  Cl 
in  the  anode  channels  is  given  by  the  conservation  equation  (5). 
According  to  the  electrochemical  reaction  (14),  the  mass  flux  of 
carbon  dioxide  is  defined  at  diffusion  layer |  membrane  interface 
as  follows: 

vSe  =  ^c°2“  (28) 

It  should  be  noted  that  Eq.  (26)  corresponds  to  an  equilibrium 
model  of  the  multi-component  mass  transfer  between  the  liquid 
and  gas  phases  in  the  anode  channels.  Derivation  of  the  auxiliary 
equation  for  the  coefficient  ifr  is  given  in  Appendix  A. 
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2.6.  Three-dimensional  model  of  DMFC 

The  three-dimensional  DMFC  model  is  based  on  the  conser¬ 
vation  equations  of  momentum,  mass  and  current.  The  coupled 
non-linear  Eqs.  (l)-(5),  (10)— (13),  (21)  and  (22)  describe  trans¬ 
fer  processes  in  the  channel,  diffusion  layer  and  membrane  on 
both  sides  of  the  fuel  cell.  Current  conservations  are  coupled 
with  diffusion  equations  via  the  interface  boundary  conditions 
that  account  for  the  electrochemical  reactions.  The  improved 
two-fluid  model  in  the  anode  channels  includes  a  new  sub-model 
for  interface  mass  transfer  Eq.  (26).  Numerical  simulation  based 
on  the  CFD  model  is  powerful  tool  for  studying  and  predicting 
the  effect  of  flow-field  geometry  on  gas  evolution,  flow  patterns 
and  fuel  cell  performance  in  DMFC  designs. 

3.  Simulation 

3.1.  Case  study  1:1.4  cm 2  DMFC  channel 

In  order  to  validate  the  three-dimensional  CFD  model  with 
the  new  sub-model  for  interface  mass  transfer,  this  case  study 
presents  simulation  results  for  a  1 .4  cm2  DMFC  channel  com¬ 
pared  with  the  conventional  sub-model  reported  by  Wang  and 
Wang  [9]  with  the  operating  conditions  given  in  Table  5. 

Fig.  1  shows  the  development  of  the  velocity  and  gas  volume 
fraction  profiles  along  the  anode  channel  at  different  cross- 
sections  from  the  inlet.  The  simulation  results  display  symmetri¬ 
cal  distributions  of  both  the  velocity  and  the  gas  volume  fraction 
in  the  channel.  The  velocity  profile  corresponds  to  a  laminar 


Table  5 

Operating  conditions  of  1 .4  cm2  DMFC 


Channel  height 

2.0  mm 

Channel  width 

2.0  mm 

Number  of  channels 

1 

Cell  length 

7  cm 

Operating  temperature 

80  °C 

Cathode  channel  pressure 

1  atm 

Anode  channel  pressure 

1  atm 

Inlet  velocity  of  anode  channel 

0.0006  ms-1 

Inlet  methanol  concentration 

1M 

Inlet  velocity  of  cathode  channel 

0.2ms-1 

Inlet  oxygen  concentration 

21  mol% 

Inlet  relative  humidity  at  cathode 

3.43  mol% 

flow  regime.  The  CFD  model  predicts  a  stable  wall  peak  of  the 
gas  volume  fraction  in  a  horizontal  channel.  The  distributions 
of  velocity  and  gas  volume  fraction  agree  with  typical  trends 
reported  in  the  literature  for  bubbly  flow  in  channels. 

The  distribution  of  methanol  and  carbon  dioxide  concen¬ 
tration  along  the  channel  are  given  in  Fig.  2.  The  methanol 
concentration  in  the  channel  is  decreased  due  to  the  electro¬ 
chemical  reaction  (14).  New  sub-model  Eq.  (26)  determines  the 
source  of  the  interface  mass  transfer  in  the  gas-liquid  flow  in 
continuity  Eq.  (4). 

Another  point  to  note  is  that  mixture  properties  in  the  dif¬ 
fusion  layer  are  functions  of  liquid  saturation.  Wang  and  Wang 
[9]  estimated  the  liquid  saturation  at  the  interface  between  the 
channel  and  the  diffusion  layer  from  an  empirical  correlation. 
The  improved  three-dimensional,  two-phase  model  eliminates 
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Fig.  1.  Distribution  of  gas  content  (£q)  and  velocity  ( u )  for  anode  channel  of  1.4  cm2  DMFC  with  operating  conditions  from  Table  1.  Average  current  density 
4500  Am“2. 
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Fig.  2.  Distribution  of  methanol  (CMe0H)  and  carbon  dioxide  concentration  (Cc°2)  for  anode  channel  of  1.4  cm2  DMFC  with  operating  conditions  from  Table  1. 
Average  current  density  4500  Am-2. 


using  an  empirical  correlation  for  liquid  saturation  (or  gas  vol¬ 
ume  fraction)  in  the  anode  diffusion  layer  (Fig.  2).  Non-uniform 
distribution  of  local  current  density  in  Fig.  3  at  the  diffusion- 
layer  |  catalyst-layer  interface  results  from  non-uniform  supply 
of  reactant  along  the  anode  channel  and  mass  transfer  in  the 
diffusion  porous  layer. 

Fig.  4  shows  the  fuel  cell  slice  with  the  predicted  distribution 
of  overpotential  ( rj  =  (Ps  —  (Pm)  at  x  =  0.03  m  from  the  channel 
inlet.  The  predictions  reveal  a  complicated  interaction  between 
the  current-collector  and  the  porous  conductor.  Strong  peaks  of 
electron  current  and  gradient  of  potential  occur  near  the  edges  of 
the  current-collector  on  both  sides  of  the  fuel  cell.  These  edges 


collect  current  under  anode  and  cathode  channels,  respectively. 
The  peak  current  is  about  10  times  higher  than  the  mean  current 
density.  Similar  results  for  a  porous  diffusion  layer  were  reported 
by  Kulikovsky  et  al.  [10]  from  a  two-dimensional  model  for  a 
gas-fed  DMFC. 

Fig.  5  provides  a  comparison  of  the  velocity,  concentration 
and  gas  content  profiles  along  the  anode  channel  calculated 
by  the  new  model  Eq.  (26)  and  profiles  corresponding  to  the 
conventional  sub-model  Eq.  (24)  with  empirical  mass  transfer 
coefficients.  The  empirical  correlation  has  been  developed  for 
mass  transfer  from  bubble  to  liquid  by  Taylor  flow  in  the  circu¬ 
lar  capillary  tube  of  a  monolithic  catalyst  reactor  [24] .  It  should 
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Fig.  3.  Distribution  of  local  electron  current  density  (/)  at  diffusion-layer  [membrane  interface  for  1.4  cm2  DMFC  with  operating  conditions  from  Table  1.  Average 
current  density  4500  A  m-2 . 
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Fig.  4.  Fuel  cell  cross-section  with  distribution  of  overpotential  (if)  at  x  =  0.03  m  from  the  inlet.  Horizontal  channel  of  1.4  cm2  DMFC  with  operating  conditions  in 
Table  1.  Average  current  density  4500  Am-2. 


be  noted  that  gas  phase  evaluation  with  mass  transfer  equations 
corresponds  to  a  non-equilibrium  model  of  transfer  processes  in 
the  anode  channels.  The  new  sub-model  Eq.  (26)  without  empir¬ 
ical  coefficients  is  in  accord  with  an  equilibrium  model  of  the 
mass  transfer  process  in  gas-liquid  flow. 

Simulation  results  in  Fig.  5  quantitatively  agree  with  the  two- 
dimensional  simulation  data  reported  by  Wang  and  Wang  [9]. 
The  void  fraction  in  the  anode  channel  increases  along  the  flow 
direction  from  0%  at  the  inlet  to  90%  at  the  outlet.  The  velocity 
increases  along  the  flow  direction  due  to  the  volume  expansion 
of  the  two-phase  mixture.  Fig.  6  presents  simulation  outcomes 
with  the  predicted  dependence  of  the  outlet  gas  content  on  the 
inlet  flow  rate  for  the  channel.  As  appears  from  the  above  three- 
dimensional  CFD  model,  an  increase  in  the  inlet  flow  rate  leads 
to  a  decrease  in  the  gas  content  at  the  outlet  section. 


Fig.  5.  Comparison  of  gas-liquid  characteristics  predicted  by  conventional  and 
new  sub-model  for  interface  mass  transfer.  Case  study  of  1.4  cm2  DMFC.  (Solid 
lines)  Averaged  three-dimensional  CFD  model  with  new  sub-model  Eq.  (26). 
(Dashed  lines)  One-dimensional  model  with  conventional  sub-model  Eq.  (24) 
reported  by  Wang  and  Wang  [9].  Average  current  density  4500  Am-2. 


Gas  content  and  channel  blocking  can  be  reduced  signifi¬ 
cantly  at  higher  liquid  flow  rates.  The  model  predictions  agree 
with  the  visualization  study  reported  by  Scott  et  al.  [1]  for  a 
transparent  DMFC  with  parallel  channels.  These  findings  are 
important  for  optimizing  the  operation  of  DMFCs. 

3.2.  Case  study  2:  DMFC  with  parallel  flow -fie  Id 

For  additional  testing  of  the  proposed  model,  simulation 
results  for  a  DMFC  are  compared  with  experimental  data 
obtained  by  Fu  and  Wang  [4]  with  the  operating  conditions  given 
in  Table  6.  They  reported  results  for  both  anode  and  cathode 
from  simultaneous  two-phase  flow  visualization  with  a  trans¬ 
parent  5  cm2  DMFC. 

Taking  operation  data  as  input,  we  have  applied  the  proposed 
model  for  the  simulation  of  a  DMFC  with  parallel  channels. 
The  distribution  of  the  calculated  velocity  and  gas  content  in  the 


Fig.  6.  Dependence  of  gas  content  on  inlet  flow  rate  of  methanol  solution.  Case 
study  of  1.4  cm2  DMFC. 


V.A.  Danilov  et  al.  /  Journal  of  Power  Sources  162  (2006)  992-1002 


999 


Table  6 

Operating  conditions  of  5  cm2  DMFC 


Channel  height 

1.5  mm 

Channel  width 

1.92  mm 

Number  of  channels 

8 

Cell  effective  area 

5.0  cm2 

Operating  temperature 

85  °C 

Cathode  channel  pressure 

15psig 

Anode  channel  pressure 

Opsig 

Inlet  methanol  stream  flow  rate 

20.8  ml  min-1 

Inlet  methanol  concentration 

2M 

Inlet  air  stream  flow  rate 

700.0  ml  min-1 

Inlet  oxygen  concentration 

21  mol% 

anode  channels  are  given  in  Fig.  7.  The  CFD  model  predicts  a 
highly  non-uniform  velocity  field  in  the  parallel  channels  and 
shows  high  values  at  the  lateral  channels  and  a  nearly  stagnant 
zone  in  the  central  region.  Simulation  results  agree  with  the 
observation  of  bubble  dynamics  in  the  channels.  The  mean  gas 
content,  sq,  along  anode  channels  computed  from  images  of 
bubbles  reported  by  Lu  and  Wang  [4]  is  between  0  and  22%. 

The  distribution  of  local  electron  current  shown  in  Fig.  8 
reveals  a  non-uniform  supply  of  reactants  through  the  diffusion- 
layer  |  catalyst-layer  interface  for  a  DMFC  with  parallel  channels. 


The  interaction  between  the  transfer  processes  in  a  DMFC 
with  parallel  flow-field  are  presented  in  Figs.  7-9.  Non-uniform 
distribution  of  flow,  gas  volume  fraction,  species  and  current 
density  profiles  in  parallel  channels  confirms  the  importance  of 
the  flow-field  in  a  DMFC  design. 

The  calculated  polarization  curve  in  Fig.  10  agrees  well  with 
the  experimental  data  for  a  5  cm2  DMFC  reported  by  Lu  and 
Wang  [4].  According  to  the  operating  conditions  in  Table  2,  the 
inlet  methanol  solution  is  saturated  with  dissolved  CO2. 

4.  Results  and  discussion 

Gas  management  depends  on  various  parameters,  e.g.,  mate¬ 
rial  properties,  cell  design  and  operating  conditions.  Channel 
blocking  restricts  the  supply  of  reactants  to  the  catalyst  layer 
and  hence  leads  to  deterioration  in  cell  electrical  performance. 
The  efficient  removal  of  carbon  dioxide  is  an  important  fac¬ 
tor  in  DMFC  design.  The  main  objective  of  gas  management 
is  to  determine  a  DMFC  design  and  operating  conditions  so  as 
to  provide  uniform  distribution  of  liquid  without  gas  accumu¬ 
lation  in  the  channels.  A  two-fluid  model  is  suitable  for  gas 
management  in  a  DMFC.  The  conventional  two-fluid  model 
uses  a  non-equilibrium  sub-model  with  empirical  coefficients 
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Fig.  7.  Distribution  of  gas  content  (sq)  and  velocity  ( u )  for  5  cm2  DMFC  with  operating  conditions  in  Table  2.  Average  current  density  2000  Am-2. 
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Fig.  8.  Distribution  of  electron  current  density  at  diffusion-layer] membrane  interface  for  5  cm2  DMFC  with  operating  conditions  in  Table  2.  Average  current  density 
2000  Am“2. 


for  interface  mass  transfer.  Wang  and  Wang  [9]  showed  the  appli¬ 
cation  of  a  homogeneous  two-fluid  model  for  a  DMFC  with  a 
conventional  sub-model  of  interface  mass  transfer  with  empir¬ 
ical  correlations.  We  suggest  a  new  equilibrium  sub-model  for 
estimating  the  interface  mass  transfer  in  anode  channels.  In  con¬ 
trast  to  Wang  and  Wang  [9],  we  use  the  improved  CFD-based 
two-phase  model  with  the  new  sub-model  for  interface  mass 
transfer  without  using  empirical  correlations.  This  study  shows 
the  application  of  the  three-dimensional  CFD  model  for  explor¬ 
ing  gas  evolution  in  a  DMFC  with  parallel  channels.  The  above 
flow-field  simulations  provide  a  good  illustration  of  the  capabil¬ 
ity  of  the  model. 

The  design  of  a  DMFC  requires  an  understanding  of  pro¬ 
cesses  such  as  mass,  momentum  transport,  electrochemical  reac¬ 
tions  and  charge  balance  that  are  taking  place  inside  the  cell. 


The  CFD  model  predicts  a  stable  wall  peak  of  the  gas  volume 
fraction  in  the  horizontal  channel.  The  calculated  distribution 
of  velocity  and  gas  volume  fraction  agrees  with  typical  trends 
reported  in  the  literature  for  bubbly  flow  in  channels.  The  cal¬ 
culated  three-dimensional  profiles  in  a  1 .4  cm2  DMFC  channel 
display  a  strongly  non-uniform  distribution  of  electronic  current 
in  the  collector  and  the  diffusion  layer.  The  current  is  concen¬ 
trated  on  the  edges  of  the  collector.  A  CFD  model  of  a  DMFC 
is  a  valuable  tool  for  developing  a  modified  collector  to  ensure 
uniform  potential  and  current  distribution. 

As  shown  in  this  study,  the  flow  geometry  of  the  anode 
side  has  an  important  impact  on  gas  evolution.  Simulation 
results  reveal  that  a  parallel  flow-field  is  not  suitable  for  gas 
management  in  a  DMFC.  Flow  maldistribution  also  occurs  in 
parallel  channels  with  single-phase  flow.  Barreras  et  al.  [25] 
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Fig.  9.  Distribution  of  methanol  concentration  at  z  =  h/2  cross-section  of  anode  channels  for  5  cm2  DMFC  with  operating  conditions  in  Table  2.  Average  current 
density  2000  Am-2. 
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Fig.  10.  Comparison  of  numerically  predicted  polarization  curve  with  experi¬ 
mental  data  for  5  cm2  DMFC.  (Points)  Experimental  data  reported  by  Lu  and 
Wang  [4].  (Solid  line)  Simulation  results  with  three-dimensional  CFD  model. 


have  performed  an  experimental  and  numerical  research  of  the 
single-phase  flow  distribution  in  a  bipolar  plate  with  a  parallel 
flow-field.  They  also  concluded  that  the  single-phase  flow  pref¬ 
erentially  moves  through  the  lateral  channels  and  results  in  an 
inappropriate  distribution  on  the  electrode  surfaces. 

It  should  be  noted  that  the  simulation  results  and  the  main  con¬ 
clusions  qualitatively  agree  with  experimental  observations  for  a 
transparent  DMFCs  available  in  literature.  One  of  the  advantages 
of  computational  modelling  is  the  ability  to  evaluate  innovative 
designs  [  1 6,26] .  The  model  developed  here  is  useful  for  the  basic 
understanding  of  three-dimensional  transport  and  electrochem¬ 
ical  phenomenon  in  DMFCs,  and  for  the  optimization  of  cell 
design  and  operating  conditions. 


5.  Conclusions 

A  CFD-based  two-phase  model  has  been  developed  for  a 
DMFC.  The  improved  two-phase  model  includes  a  new  sub¬ 
model  for  interface  transfer  without  using  any  empirical  cor¬ 
relation  for  gas-liquid  flow  in  the  anode  channels.  Simulation 
results  agree  with  typical  trends  for  bubbly  flows  in  channels. 
It  is  found  that  flow-field  design  has  a  very  significant  influ¬ 
ence  on  to  the  performance  and  gas  management  of  DMFCs. 
The  CFD  model  is  validated  against  the  polarization  curve  for  a 
5  cm* 1 2 3 4 5 6 7 8 9  DMFC.  An  increase  in  inlet  methanol  solution  flow  rate 
is  found  to  lead  to  a  decrease  in  the  total  outlet  gas  content  on 
the  anode  side.  Interaction  between  the  current  collector  and 
the  diffusion  layer  leads  to  a  complex  distribution  of  current  and 
potential  in  fuel  cell.  Computational  results  furnish  explanations 
of  the  observed  flow  in  transparent  DMFCs,  as  reported  by  other 
workers.  The  proposed  model  is  valuable  in  CFD-based  DMFC 
design. 
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Appendix  A 


Assuming  the  coefficient  xj/  is  constant,  we  can  evaluate  it 
from  total  mass  balance  for  gas  phase 

6rin  —  Gout  =  TCjVanode  (A.l) 


where  Gm  is  the  gas  flow  rate  at  inlet  section;  Gm  =  0;  Gout  the 
gas  flow  rate  at  the  outlet  section;  Vanode  is  the  volume  of  anode 
channels.  The  gas  flow  rate  at  the  outlet  section  is 

Gout  =  ^G,outPG^G  ^out  (A. 2) 


where  Sout  is  the  area  of  the  stream  outlet  section;  uc,out  the 
velocity  of  gas  at  section  Sout;  £°ut  is  the  gas  content  in  the  two- 
phase  stream  at  section  Sout.  The  gas  content  at  the  outlet  section 
of  anode  channel  is  calculated  using: 

gOUt  _  _ K°UVmol,L _  ,A 

°  K°UtAnol,L  +  (1  -  Kout)Pmol,G 

The  mean  source  of  gas  is: 

r0  =  M0j>^  +  ^)  (A.4) 

where  Mq  is  the  molecular  weight  of  gas  phase;  xj/  the  coef¬ 
ficient,  x[/  =  L/bV ;  y  the  mean  fractional  vaporization,  y  = 
(1/ Vanode)  f  V  dV.  The  mean  mass  flux  of  the  carbon  dioxide 
component  in  the  anode  channel  is: 

=  (A.5) 

The  fractional  vaporization  yout  at  the  outlet  section  of  the  anode 
channel  is  given  by: 


X(Kj  -  1  )C°ut(Mmix/M,-) 

^  (K,  -  l)yout  +  1 


where  component  molar  fraction  at  the  outlet  section  is  evaluated 
from  mass  balance. 

Finally,  solving  Eqs.  (A.4)  and  (A.l)  for  x[/  we  have  [16,26]: 


xj/  = 


Gin  —  G 


out 


Tanod 


McoUi 


(A.7) 
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